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Abstract: Motivated by ongoing experimental analyses, we report on the calculation of 
Ci^ next-to-leading order QCD corrections to the production of tt pairs in association with two 

p I hard b-jets at the Fermilab TeVatron. Besides the total cross section and its scale depen- 

dence, a few differential distributions applicable for Higgs boson searches at the TeVaron 
are given. The QCD corrections with respect to leading order are negative and small. For 
our main setup they amount to 2%, and remain reasonably stable against changes of cuts. 
This proves that an integrated next-to-leading order /C-factor does not necessarily need to 
be applied in the background estimation for the ttH — >• ttbb signal process. The distribu- 
tions show similarly small corrections. The shape of kinematic distributions is distorted at 
most by about 20% in some corners of the phase space. Even though it is not the main 
purpose of this paper, we also evaluated the forward-backward asymmetry of the top quark 
at next-to- leading order. 
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1 Introduction 

The production of the Standard Model Higgs boson in association with a top-anti-top 
pair allows for a direct study of the top Yukawa coupling, which is an important step 
in understanding the nature of the electroweak symmetry breaking mechanism. The tiH 
production channel is especially important in the range of Higgs boson masses mn < 
140 GeV, where the Higgs boson decays predominantly into bb pairs. While the LHC 
is making ground breaking progress towards discovering the Higgs boson, and certainly 
tightens exclusion bounds well past those obtained at the TeVatron, the analysis of the data 
from the latter is still being completed. It is a fact that the predicted cross section for Higgs 
boson radiation off top quarks at the TeVatron is rather low [1-4]. Therefore, a discovery of 
the Higgs boson in this channel alone would not have been possible. However, it contributes 
to the combination of the Standard Model Higgs boson searches [5]. The main background 
process consists of direct production of the final state without resonances, i.e. of the QCD 
generated process ttbb. The procedure used by the experimental collaborations for the 
background estimate, is either direct use of leading order (LO) Monte Carlo simulations 
[6], or additional reweigthing by a K-factor [7]. Unfortunately, a next-to- leading order 
(NLO) K-factor for this process under TeVatron conditions has never been published. 
According to [7], the analysis employs the published value [8] for the LHC. Since the 
production mechanisms are very different in both cases, it is questionable that such a 
procedure leads to reliable estimates. In this paper, we address this issue by providing an 
NLO QCD prediction to ttbb, the irreducible QCD background to the ttH — t- ttbb process 
at the TeVatron. 
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2 Theoretical Framework 



At tree level, tt production in association with two b-jets proceeds via the scattering of two 
gluons or two quarks. A few examples of LO graphs contributing to pp — t- ttbb production 
are shown in Figure 1. The virtual corrections are obtained from the interference of the sum 
of all one-loop diagrams with the Born amplitude. One can classify them into self-energy, 
vertex, box-type, pentagon-type and hexagon-type corrections. In Figure 2 a few examples 
of pentagon and hexagon diagrams contributing to the virtual corrections to the pp — )• ttbb 
process are given. And finally, the real emission corrections to the LO process arise from 
tree level amplitudes with one additional parton, an additional gluon, or a quark anti-quark 
pair replacing a gluon. All possible contributions can be divided into four subprocesses, 
qq — )• ttbbg, gg — )• ttbbg, qg — )■ ttbbq and gq — )■ tibbq. In Figure 3 a representative set of 
Feynman diagrams contributing to the real emission corrections is shown. 

The calculation of NLO corrections to pp — t- ttbb production proceeds along the same 
lines as our earlier work on pp — t- ttbb [8], pp{pp) — )• tijj [9, 10] and pp{pp) — ?• W^W~bb — )• 
e'^Ve^Ji~y^ibb [11]. The methods and internal tests developed there have therefore been 
straightforwardly adapted for this project. Let us stress here, that results for the pp — )• ttbb 
and pp{pp) — )• W~^W~bb processes have also been obtained by other groups [12-14]. 

To summarize briefly, off-shell methods and the OPP reduction procedure [15], as 
implemented in the Helac-NLO system [16] are used in computing the NLO QCD cor- 
rections for the pp —7- ttbb processes. The system consists of, CutTools [17-20] and 
Helac-ILoop [21], which handle the virtual corrections and Helac-Dipoles [22] for the 
real emission contributions. For the phase space integration, the Kaleu package [23] is 
used and results are cross checked with the help of the Phegas phase space generator [24]. 
Moreover, OneLOop [25] is employed for the evaluation of the one-loop scalar functions. 

LO results, which are generated with Helac-Dipoles have been cross checked with 
Helac-Phegas^ [26-28]. Perfect agreement has been found. The numerical stability of 
virtual corrections has been monitored by checking Ward identities at every phase space 
point. Those events, which violate gauge invariance (0.13% out of all 2 X 10^ generated 
unweighted events for which one-loop contributions are calculated) have been recalculated 
with quadruple precision. In addition, the cancellation of soft and collinear divergences 
after combining virtual and real corrections has been checked numerically for a few phase 
space points. 

3 Phenomenological Results 

In the following, we present predictions for the tibb+X process at the TeVatron run II with 
^/s = 1.96 TeV. We use the Tevatron average mass of the top quark mt = 173.3 GeV as 
measured by the CDF and DO experiments [34]. The masses of all other quarks, including 
b quarks, are neglected. We have consistently employed the MSTW2008 set of parton 
distribution functions (PDFs) [35]. In particular, we take MSTW2008LO PDFs with 1-loop 
running at LO and MSTW2008NLO PDFs with 2-loop running Us at NLO, including 

^Let us stress that Helac-Phegas has already been extensively used and tested, see e.g. [29-33]. 
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Figure 1. A representative set of Feynman diagrams contributing to the leading order pp — > ttbb 
process. Double lines correspond to top quarks, single lines to light quarks and wiggly ones to gluons. 




Figure 2. A representative set of pentagon and hexagon diagrams contributing to the virtual cor- 
rections to the pp ttbb process. Double lines correspond to top quarks, single lines to light quarks 
and wiggly ones to gluons. 




Figure 3. A representative set of Feynman diagrams contributing to the real emission corrections 
to the pp ttbb process. Double lines correspond to top quarks, single lines to light quarks and 
wiggly ones to gluons. 
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five active flavors. Contributions induced by bottom-quark densities are not taken into 
account due to their negligible size. The renormalization and factorization scales are set to 
a common value /i/j = fip = = mt- All final-state massless partons with pseudorapidity 
|?7| < 5, defined as rj = — In [tan(^/2)], where 9 is the angle between the parton momentum 
and the beam axis, are recombined into jets with a resolution parameter R = 0.8 via an 
IR-safe algorithm. We have applied three different jets algorithms: [36-38], anti-kT 
[39] and the Cambridge/ Aachen (C/A) algorithm [40]. For our main setup, two b-jets are 
required to have 

pxib) > 20 GeV, \y{b)\ < 2.5, Ai?^^ > 0.8 , (3.1) 

where PT{b), y{b) are the transverse momentum and rapidity of the b-jet, and Aii^j is the 
separation in the plane of rapidity and azimuthal angle between bb pairs. Jets momenta 
are formed as the four- vector sum of massless parton momenta. At LO, there are exactly 
two massless final state partons, which are identified as two b-jets, provided they pass 
the cuts described above. At NLO, a third parton might emerge. It could be recombined 
with another parton to give a b-jet, or an additional jet with unrestricted kinematics may 
appear. Outgoing top and anti-top quarks are left on-shell, they do not undergo any cut 
selection. 

3.1 Integrated Cross Sections 

We start with the total cross sections. In Table 1 and Table 2 integrated cross sections 
at LO and NLO for pp — )• ttbb + X production at the TeVatron run II are presented. In 
Table 2, the scale dependence of the total cross section is also given. Using a fixed scale 
independent of the final state kinematics, which we set to mt and estimate the error with 
the usual variation in the range between nit/ 2 and 2mt our findings can be summarized as 
follows 

(TLo (TeVatron, mj = 173.3 GeV, MSTW2008lo ) = 3.912 +J;^o5jJ35j fb , (3.2) 

f^NLO (TeVatron, mt = 173.3 GeV, MSTW2008nlo) = 3.835 ^? 015(26%) ^ ' ('^••^) 

which leaves us with an NLO /C-factor equal to /C = 0.98 and a negative NLO QCD 
correction of the order of 2%. This is very different from the LHC case, where the same 
cut selection implies NLO corrections of the order of 77% [8, 12]. There, a dynamical scale, 
/i = 'mt\/pT{b)pT(b), and m^j > 100 GeV had to be introduced in order to reduce the 
high corrections down to 25% — 30% [13]. The dissimilarity between the size of NLO QCD 
corrections comes mostly from the difference in the production process. At the TeVatron, 
with our cut selection the qq channel dominates the total LO pp cross section at about 91% 
followed by the gg channel with about 9%. In contrast, the gg channel comprises about 
94% of the LO pp cross section at the LHC, followed by the qq channel with about 6%. 

To assess the effect of changing the jet algorithm, we compare NLO results for different 
jet finders, the k^, anti-kT and the inclusive Cambridge/ Aachen (C/A) jet algorithms as 
shown in Table 1. No significant change of the results due to the choice of the jet algorithm 
is observed. Differences are below 1%. 
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In Equation (3.2) and Equation (3.3) the scale dependence is indicated by the upper 
and lower values. The upper (lower) value represents the change when the scale is shifted 
towards ^ = mt/2 (/U = 2mt). Rescaling the common scale from the default value mj up 
and down by a factor 2 changes the cross section at LO by 89%. On the other hand, the 
improvement in the scale stability at NLO is prominent. The scale uncertainty is reduced 
down to 26%. 





CTLO [fb] 


anti— kT rrui 


4lo [fb] 


^NLO [fb] 


Q^max — 1 


3.912(3) 


3.835(3) 


3.859(3) 


3.853(3) 


"max = 0.01 


3.912(3) 


3.836(5) 


3.861(5) 


3.856(5) 



Table 1. Integrated cross section at LO and NLO for pp — > ttbb + X production at the TeVatron 
run IL Results for three different jet algorithms and two different values of Ofmax presented. The 
scale choice is ^ = rut- 





0.5 • iJit 


1 • mt 


2 • mt 


CTLO [fb] 


7.408(5) 


3.912(3) 


2.207(2) 


0-NLO [fb] 


4.827(8) 


3.835(3) 


2.820(3) 



Table 2. Scale dependence of the total cross section for pp — > ttbb -\- X production at the TeVatron 
run II at LO and NLO with fi ^ ^ ■ mt- 

In addition, integrated NLO cross sections for two values of the unphysical cutoff 
parameter amax) which is a common modification of subtraction terms in the phase space 
away from the singularity first introduced in [41], are given in Table 1. To be specific, 
Omax = li which Corresponds to the original formulation of [42, 43], and Omax = 0.01 are 
considered. The independence of the final result on the value of the amax parameter is 
obtained at the per-mil level. This is a strong consistency check of the calculation of the 
real emission part. For more details on the cimax implementation in the Helac-Dipoles 
package see e.g. [8, 22]. 

Results for a slightly modified setup have also been generated, to determine the stability 
of the NLO /C-factor. On the one hand, a higher transverse momentum cut on the b-jet 
of 40 GeV has been chosen. On the other hand, a higher jet separation cut, ARf^i > 1 has 
been used together with a new jet resolution parameter R = 1. All other parameters have 
been left unchanged. With this modified selection of cuts the integrated NLO /C-factor has 
changed from 0.98 to 0.88. More precisely, we have obtained the following integrated cross 
section at LO and NLO 

o-Lo = 0.8135(6) fb , (3.4) 

CTNLO = 0.7121(10) fb, (3.5) 
which results in negative NLO QCD corrections of the order of 12%. 
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Figure 4. Distribution of the invariant mass m^;, (left panel) and distribution in the transverse 
momentum pt^^^ (right panel) of the bottom-anti-bottom pair for pp ^ tibb + X at the TeVatron run 
II at LO (blue dashed line) and NLO (red solid line). The lower panels display the differential K, 
factor. 




Figure 5. Distribution in the rapidity y^^ (left panel) of the bottom- anti-bottom pair and distribution 
of the ARjjI separation (right panel) for pp — >■ tibb -\~ X at the TeVatron run II at LO (blue dashed 
line) and NLO (red solid line). The lower panels display the differential K, factor. 

3.2 Differential Cross Sections 

In the following, the impact of QCD corrections on differential cross sections is analyzed. 
The differential distributions relevant for Higgs boson searches in the tiH — t- tibh channel 
are plotted first. In Figure 4, the distribution of the invariant mass m^i and the distribution 
in the transverse momentum pT^- of the bottom-anti-bottom pair for pp — t- tibb -|- X at the 
TeVatron run II is plotted. The dashed curve corresponds to the LO, whereas the solid one 



-6- 




Figure 6. Distribution in the transverse momentum px^ (left panel) and distribution in the rapidity 
Ub (right panel) of the bottom quark for pp — > ttbb + X at the TeVatron run II at LO (blue dashed 
line) and NLO (red solid line). The lower panels display the differential K, factor. 

to the NLO result. The upper panels show the distributions themselves while the lower 
panels display the ratio of the NLO value to the LO result, calculated according to 

for an observable O under investigation, called the differential or dynamical KL factor. In 
Figure 5, the distribution in the rapidity of the bottom-anti-bottom pair and the distri- 
bution of the Ai2j5 separation between two b-jets are shown. In all cases, the small NLO 
corrections of the integrated cross section are also visible at the differential level. In case of 
the most important observable, the m^^i differential distribution, in the phenomenologically 
relevant region, i.e. below 140 GeV, corrections of the order of ~ 10% are reached. Besides, 
the separation between the b-jets of 0.8 together with the transverse momentum cut on 
b-jets of 20 GeV sets an effective lower bound on the invariant mass of two b-jets of the 
order of m^j > 15.6 GeV [8]. If we had chosen Ai?^^ > 0.5 the minimum m^^ would have 
been around 9.9 GeV instead. 

In view of the small NLO corrections to the total cross sections as well as to the 
differential distributions presented here, we conclude that a meaningful analysis at the 
TeVatron can be performed with the present setup, i.e. with the renormalization and 
factorization scales fixed to a common value (the mass of top-quark). 

Finally, for completeness the distributions for b-jet kinematics are given in Figure 6. 
Namely, the distribution in the transverse momentum px^ and the distribution in the rapid- 
ity of the bottom quark are shown. Similarly to the bottom-anti-bottom pair kinematics, 
also here, small NLO corrections are visible. The shape of kinematic distributions can be 
distorted by 20% at most in some regions of phase space. This is once again contrary to 
the LHC case, where large and relatively constant NLO corrections have been obtained for 
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-1.5 -1 -0.5 0.5 1 1.5 -1.5 -1 -0.5 0.5 1 1.5 

yt yt 

Figure 7. Differential cross section distributions as a function of rapidity, yt, of the top and anti- 
top quark at LO (left panel) and NLO (right panel) for pp — >■ ttbb + X production at the TeVatron 
run II. The (orange) solid curve corresponds to the top quark, whereas the (brown) dash-dotted one 
to the anti-top quark. 
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yt yt 



Figure 8. Differential asymmetry A{t), as a function of the top quark rapidity at LO (left panel) 
and NLO (right panel) for pp — > ttbb + X production at the TeVatron run II. NLOi refers to a 
result with a consistent expansion in as , while NLO2 to the unexpanded one. 

fi = mt, whereas a dynamical scale and additional cuts allowed to reduce them down to 
20% - 40%. 

3.3 Forward-Backward Asymmetry 

As a bonus of our study, we calculated the integrated top quark forward-backward asym- 
metry for the ttbb production process at the TeVatron. At LO the asymmetry is defined 
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as 

.t _ cTLoiyt > 0) - (TLojyt < 0) „^ 
^^^'"°"cTLo(y*>0)+aLo(yt<0)' ^'-'^ 

where yt is the rapidity of the top quark and 0-^0 = (^Loiut > 0) ibo"Lo(?/i < 0) is evaluated 
with LO PDFs and LO a^. On the other hand, the asymmetry at NLO is expressed through 

At _ '^LO + '^^NLO (o o\ 

where Sa^j^Q are the NLO contributions to the cross sections and are evaluated this 
time with NLO PDFs and NLO a^. The ratio generates contributions of 0{a1) and higher, 
which are affected by unknown next-to-next-to-leading order contributions. Therefore, it 
is necessary to expand Equation (3.8) to first order in a^. The following definition is then 
obtained [44, 45] 

At _ ^LO ( 1 j_ '^'^NLO ^^NLO \ (n Q\ 

The integrated forward-backward asymmetry of the top quark at LO for tibh produc- 
tion amounts to (the error in parentheses corresponds to scale variation) 

^FB,LO = -0-088(2). (3.10) 

With NLO corrections, the asymmetry is reduced down to 

-^FB.NLO = -0.044(6) , (3.11) 

when a definition with a consistent expansion in as is used. For an unexpanded ratio of 
the NLO cross sections, the result is 

^FB.NLO = -0.061(14) . (3.12) 

The two NLO definitions give results which differ by about 40% for the central scale. Also 
the theoretical error as calculated from the scale dependence is more than a factor 2 higher 
in the latter case. 

In Figure 7, the rapidity distributions for the top and anti-top quarks are presented 
at LO and NLO. Results are not symmetric around yt = and are shifted to a forward 
direction for the anti-top quarks and a backward direction for the top quarks. This shows 
that anti-top quarks are preferentially emitted in the direction of the incoming protons. 

In Figure 8, we have also plotted the differential asymmetry, A{yt)- It rises up to 
±20% at LO and well above ±10% at NLO in suitably chosen kinematical regions. After 
including the NLO corrections the forward-backward asymmetry of top quarks is reduced 
by a factor 2 as can be seen both from Equation (3.10) and Equation (3.11) as well as in 
Figure 8. 
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4 Summary and Conclusions 

In this paper, a computation of the NLO QCD corrections to the top quark pair production 
in association with two hard b-jets at the TeVatron run II has been presented. The total 
cross section and its scale dependence, together with a few differential distributions have 
been given. The impact of the NLO QCD corrections on the integrated cross sections as 
calculated for the fixed scale jip = tJ-R = M = iTT-t, is small, of the order of 2%. More- 
over, the NLO QCD corrections reduce the scale uncertainty of the total cross sections as 
compared to LO calculations. As a further matter, the NLO corrections to the differential 
distributions are below 10% in the phenomenologically significant regions, reaching 20% at 
some corners of the phase space. 

Since at the TeVatron, the corrections are small, the integrated and differential /C- 
factors do not necessarily need to be applied in the background estimation for the ttH — >• 
ttbb signal process at the TeVatron. At least not for the observables which have been scru- 
tinized here. Definitely, an application of the large LHC /C-factor to TeVatron phenomeno- 
logical analyzes dramatically overestimates the size of the irreducible tibb background. 

Let us conclude by noting that the predicted rates are very small. Combined with 
the luminosity of the TeVatron they amount to just a few events. The relevance of the 
present study lies, however, in the fact that Higgs exclusion bounds are obtained using a 
combination of data from different production channels with the ttH process, and thus its 
irreducible background, also taken into account by the CDF and DO collaborations. 
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